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Abstract 

Background: The Anaphase Promoting Complex or Cyclosome (APC/C) is the largest member of the ubiquitin 
ligase [E3] family. It plays a crucial role in the control of the cell cycle and cell proliferation by mediating the 
proteolysis of key components by the proteasome. APC/C is made of a dozen subunits that assemble into a large 
complex of -1.5 MDa, which interacts with various cofactors and targets. 

Results: Using comparative genomic and phylogenetic approaches, we showed that 24 out of 37 known APC/C 
subunits, adaptors/co-activators and main targets, were already present in the Last Eukaryotic Common Ancestor 
(LECA) and were well conserved to a few exceptions in all present-day eukaryotic lineages. The phylogenetic 
analysis of the 24 components inferred to be present in LECA showed that they contain a reliable phylogenetic 
signal to reconstruct the phylogeny of the domain Eucarya. 

Conclusions: Taken together our analyses indicated that LECA had a complex and highly controlled modern-like 
cell cycle. Moreover, we showed that, despite what is generally assumed, proteins involved in housekeeping 
cellular functions may be a good complement to informational genes to study the phylogeny of eukaryotes. 
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Background 

The anaphase-promoting complex or cyclosome (APC/ 
C) has been recently characterized as a member of the 
ubiquitin ligase family (also referred as E3) (for a recent 
review see [1]). E3s mediate the transfer of one or sev- 
eral ubiquitin monomers on a protein substrate in a 
two-step reaction involving at least three partners. First, 
an ubiquitin-activating enzyme (El) activates and trans- 
fers ubiquitin to an ubiquitin-conjugating enzyme (E2). 
Next, E3 mediates the transfer of ubiquitin from E2 to a 
lysine residue of the target protein. Both steps require 
ATP. Most E3s are able to polyubiquitinate proteins by 
adding new ubiquitin monomers to the first attached 
one [2]. Polyubiquitinated proteins are targeted to the 
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26S proteasome for degradation, whereas mono-ubiqui- 
tinated proteins can be altered in their function or sub- 
cellular location by proteins containing ubiquitin- 
binding domains [2]. E3s are divided in several families 
according to the presence of signature motifs. Among 
them, E3s containing a HECT (Homologous to E6-AP 
Carboxyl Terminus) domain receive ubiquitin from E2 
before attaching it on the substrate, whereas E3s harbor- 
ing a RING (Really Interesting New Gene) finger- 
domain mediate the transfer of ubiquitin directly from 
E2 to the substrate [2,3]. RING finger E3s form the lar- 
gest family and may also contain a subunit with a Cullin 
domain [1]. Among them, the APC/C is atypically large 
(1.5 MDa) and complex, being composed of one or sev- 
eral copies of at least a dozen subunits and of various 
adaptors/co-activators (Figure 1) [4]. The function of 
the APC/C has been extensively studied in animals and 
yeast, where it was shown to have a critical role in cell 
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Structural complex 

Figure 1 Schematic structure of the APC/C. The APC/C is divided in three parts: the structural complex composed of three subunits (in pink), 
the TPR arm made of nine proteins (in blue) and a second arm (the catalytic arm) made of three proteins (in green) interacting with the E2 
protein (in purple) and the substrate (in grey). The activity and specificity of the APC/C is modulated by various adaptor/co-activators (in orange). 
The location of Apc14 (in black) remains uncertain. The presence of RING finger and Cullin conserved functional domains, and of TPR and WD 
repeats is indicated. 



cycle progression through the tight control of degrada- 
tion of key proteins (e.g. mitotic cyclins, anaphase inhi- 
bitors, spindle proteins, regulatory kinases, inhibitors of 
DNA replication, etc.) [5,6]. 

Electron microscopy observations, in vitro assays, 
genetic experiments and structural studies have shed 
light on the composition, structural organization, assem- 
bly and molecular activity of the APC/C (Table 1 and 
Figure 1) [1,4,7]. The APC/C core is divided in three 
functional parts: i) the structural complex, which is 
made of Apcl, Apc4 and Apc5 subunits (in pink, Figure 
1), serves as scaffold; ii) the catalytic arm (in green, Fig- 
ure 1) that houses the E2-binding site is made of ApclO 
and of the Apc2 and Apcll proteins that contain the 
Cullin and RING finger domains, respectively; and iii) 
the TPR arm (in blue, Figure 1) allows positioning both 
the E2 and the substrate in order to promote the ubi- 
quitin transfer. This second arm is composed of four 
subunits (Apc3, Apc6, Apc7, and Apc8) containing tet- 
ratricopeptide repeats (TPR) and of several accessory 



proteins (Apc9, Apcl2, Apcl3, and Apcl5, the last two 
being essential in Schizosaccharomyces pombe [8]). It 
has been proposed that the TPR repeats interact with 
the isoleucine- and arginine-rich motifs found in the C- 
terminal regions of adaptors/co-activators (see [1] and 
references therein for additional details). The TPR arm 
may also contain Apcl6, a subunit recently reported by 
the MitoCheck consortium [9]. Finally, the location of 
Apcl4, a yeast essential subunit [8] remains 
undetermined. 

The APC/C activity and specificity are modulated by 
several adaptors/co-activators (Cdc20, Cdhl, Amal, 
Mfrl, Rap and Cortex) (Table 1 and Figure 1) [7,10]. 
These are paralogous proteins containing WD-repeats 
that mediate the interaction between the APC/C and 
the D-, KEN-, A- or O-boxes present on target sub- 
strates (see [11] and references therein). Among those 
adaptors, Cdc20 and Cdhl are the most important, 
being directly involved in the activation and substrate 
selectivity of the APC/C at different stages of the cell 
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cycle (reviewed in [1,6,7]). The interaction of the APC/C 
and either Cdc20 (in early-mid-mitosis) or Cdhl (during 
late-mitosis and Gl-S transition) is strongly dependent 
on the high or low activity of Cdks (Cyclin dependent 
kinases) [6]. Briefly, Cdc20 activates the APC/C during 
early mitosis once the chromosomes are properly 
attached and bi-oriented at the metaphase plate during 
a process known as the spindle assembly checkpoint 
[12]. The APC/C cdc20 targets securins and cyclins Bl 
towards destruction by the proteasome. The degradation 
of these two proteins promotes the activation of 
separases, which then cleave the cohesin complex (CC) 
leading to the separation of sister chromatids and the 
initiation of the anaphase [6,12,13]. During anaphase, 
the APC/C cdhl targets Polo-like kinase 1, Aurora 
kinases, mitotic cyclins and Cdc20 towards degradation 
leading to the exit of mitosis. The APC/C" cdhl remains 
active during the Gl/S phase ensuring the degradation 
by the 26S proteasome of several inhibitors of DNA 
replication, thus allowing the synthesis of DNA [7]. At 
the end of the S phase, the increase of the activity of 
Cdks inhibits the interaction between Cdhl and the 
APC/C complex, precluding new rounds of DNA synth- 
esis [1]. By contrast, other APC/C activators seem to 
have more restricted roles: Amal is required for sporu- 
lation and during the anaphase of meiosis I in budding 
yeast [7,14]; Mfrl acts at the end of meiosis II in S. 
pombe [10]; Cortex encodes a putative Drosophila mela- 
nogaster female meiosis-specific co-activator of the 
APC/C prior to the metaphase I arrest [15] and, finally, 
Rap (retina aberrant in pattern) mediates the degrada- 
tion of cyclins during the development of eye imaginal 
discs in D. melanogaster [16]. 

If most of the APC/C studies have been carried out in 
yeast and animals, recent experiments with the land 
plant Arabidopsis thaliana have allowed the identifica- 
tion of 12 transcribed genes that are homologous to ver- 
tebrate and yeast APC/C subunits and of eight Cdhl/ 
Cdc20 homologues (Table 1) [17-20]. By contrast, very 
little information is available for representatives of the 
other major eukaryotic lineages. The only exception 
concerns the kinetoplastid species Trypanosoma brucei, 
shown to encode seven APC/C subunit homologues in 
its genome (Table 1) [21]. The apparent conservation of 
components of the APC/C in these few distantly related 
eukaryotes opens the question of the origin and evolu- 
tion of this atypically large and complex ubiquitin-ligase. 
To tackle this issue, we carried out an in-depth analysis 
of 16 APC/C subunits and six adaptors/co-activators 
(Table 1) in all eukaryotic lineages for which representa- 
tives with complete (or nearly complete) genome 
sequences were available (Figure 2). We also included in 
our study several major direct or indirect targets of 
APC/C, namely the separase, the securin, cyclins A and 



B, Cdks-1 and -2 and the nine components of the cohe- 
sin complex (Table 2). The phylogenomic analysis of 
these proteins supports that most of them were present 
in the last eukaryotic common ancestor (LECA), indicat- 
ing that this organism likely possessed a highly con- 
trolled cell cycle that may have been very similar to that 
of present-day eukaryotes. Finally our analyses indicate 
that APC/C components and targets carry a bona fide 
phylogenetic signal that can be used to trace back the 
evolutionary history of the eukaryotic domain. 

Results and Discussion 

Most APC/C components and main targets were present 
in LECA 

We used a phylogenomic approach [22] in order to 
study the origin and evolution of APC/C and its main 
targets in eukaryotes. The first step consisted in the sur- 
vey of complete (or nearly complete) genome sequences 
available in public databases to retrieve homologues of 
each component of this system. Working on complete 
genomes ensures the rigorous inference of the presence 
or absence of homologues in each genome. Then, phylo- 
genetic analyses allow inferring the origin and the subse- 
quent evolution of each component. 

We searched for orthologues in 65 taxa representing 
the eukaryotic diversity (all major groups were present 
except Rhizaria and Cryptophyta). More precisely, our 
taxonomic sampling covered (i) Holozoa (represented by 
Metazoa, and their close relatives the Choanoflagellata 
and Capsaspora), (ii) Fungi (with three main lineages: 
Dikarya, Chytridiomycota and Microsporidia) and (iii) 
Apusozoa. Whereas the position of Apusozoa remains 
uncertain, Metazoa, their unicellular allies and Fungi 
represent the opisthokont lineage that together with 
Amoebozoa, form one of the two putative major divi- 
sions of eukaryotes, the Unikonta [23]. The other major 
division, the Bikonta, was represented in our study by 
genomes from Excavata, Heterokonta, Apicomplexa and 
Ciliata (two main groups within the Alveolata), Hapto- 
phyta, and Viridiplantae and Rhodophyta (two main 
Plantae lineages) (Figure 2). 

At this step, it is interesting to notice that, except for 
adaptors/co-activators and a few other exceptions, we 
identified at most only one homologue of each APC/C 
component and main target coding genes in each gen- 
ome (Additional file 1, Table SI and Additional file 2, 
Table S2). In addition, some of them were found only in 
very restricted sets of species. For example, orthologues 
of Apcl4 and of two subunits of the TPR arm (Apc9 
and Apcl5) were present only in a few ascomycetes 
(Table 1 and Additional file 1, Table SI) suggesting that 
they are recent innovations that emerged after the first 
diversification of fungi (Figure 2). The TPR arm protein 
Apcl6 may be even more recent because it was present 
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Figure 2 Evolutionary history of APC/C subunits, adaptors/co-activators and main targets We used a reference phylogenetic tree rooted 
in-between Unikonta and Bikonta [23] showing the relationships between the eleven major eukaryotic lineages for which complete (or nearly 
complete) genome sequences are available [56,63,64,78]. Holozoa are representated by Capsaspora, Choanoflagellata and Metazoa (in pink); 
Fungi by Dikarya, Chytridiomycota and Microsporidia (in dark purple); Alveolata by Apicomplexa and Ciliata (in brown) and Plantae by 
Viridiplantae and Rhodophyta (in green). The completion status for each genome sequence is indicated: "C": complete, " A": draft assembly, "P": 
ongoing (accordingly the loss of some components inferred in the corresponding lineages should be confirmed when complete genome 
sequences are available). LECA (grey circle) represents the Last Eukaryotic Common Ancestor. Based on phylogenetic analyses, we inferred the 
time of appearance or loss of orthologues of APC/C subunits (dark blue), activators (light blue) and main targets (pink). For each lineage, protein 
gains and losses are indicated by green and red arrows, respectively. 
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Table 2 Taxonomic distribution of homologues of major APC/C targets in major lineages of eukaryotes. 

Core complex Loading complex Cohesion establishment Maintenance Dissolution 
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>50 
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100, > 50, and < 50 indicate that at least one orthologue has been detected respectively in all, at least 50%, and less than 50% of the studied genomes for a 
given lineage. 



only in the two Gnathostomata representatives within 
metazoa (Table 1 and Figure 2). By contrast, based on 
ML and Bayesian phylogenetic analyses (Additional file 
3, Data SI), we identified orthologues of the other 12 
APC/C subunits and of two adaptors/co-activators in at 
least two bikont and two unikont major lineages, indi- 
cating that they were likely present in LECA (Table 1 
and Figure 2). Accordingly, their absence in any pre- 
sent-day genome should be interpreted as secondary 
loss. This is for instance the case of the six subunits 
composing the catalytic arm and the structural complex, 
for which we detected homologues in all opisthokonts 
but choanoflagellates, which ancestrally missed Ape 11 
and Apc5, suggesting ancient losses in this lineage. Simi- 
larly, whereas the four main subunits of the TPR arm 
(Apc3, Apc6, Apc7 and Apc8) were inferred to be pre- 
sent in the ancestor of opisthokonts, Apc7 was missing 
in ascomycete and basidiomycete representatives sug- 
gesting a secondary loss in the ancestor of these two 
fungal groups after their separation from chytrids. 
Regarding the other proteins associated to the TPR arm, 
beside the case of Apc9, Apcl5 and Apcl6 already men- 
tioned, homologues of Apcl2 and Apcl3 were poorly 
represented in opisthokonts. More precisely they were 
missing in choanoflagellates, Capsaspora, most fungi 
and some animals. However, the presence of Apcl2 and 
Apcl3 orthologues also in some bikont lineages indi- 
cated that these subunits were present in LECA and 
thus in the ancestor of opisthokonts. Accordingly, their 
poor taxonomic distribution in this eukaryotic lineage 
results from convergent secondary losses (Figure 2). 
Finally, orthologues of the two adaptors Cdc20 and 



Cdhl were present in all opisthokonts. Among them, 
the case of Microsporidia deserved attention. Indeed, 
whereas only one APC/C subunit (ApclO, located in the 
catalytic arm) had previously been reported in Encepha- 
litozoon cuniculi [24], we additionally found orthologues 
of one component of the structural complex (Apcl), 
three of the TPR arm (Apc3, Apc6 and Apc8) and of 
two adaptors/co-activators (Cdhl and Cdc20) in gen- 
omes of four representatives of this group of highly 
derived parasitic anaerobic fungi. The conservation of at 
least one component of each functional part of the 
APC/C suggested that a minimalist version of the APC/ 
C might exist in Microsporidia. More drastic losses 
were observed in the anaerobic parasite Entamoeba his- 
tolytica where the absence of all but four components 
contrasted with the conservation of all 14 subunits and 
adaptors/co-activators inferred to be present in LECA in 
the second amoebozoan studied (the free-living Dictyos- 
telium discoideum) (Table 1 and Figure 2). Such massive 
losses were also observed for the parasitic excavate 
Giardia intestinalis. However, in contrast with Micro- 
sporidia, the more reduced set of components and, 
more precisely, the absence of all proteins composing 
the TPR arm appeared less compatible with a fully 
operational APC/C system in these two anaerobic 
parasites. 

In bikonts, orthologues of the 12 components inferred 
to be present in LECA were also inferred to be present 
in the ancestors of Plantae and Heterokonta (Figure 2). 
However, in red algae, the haptophyte Emiliania huxleyi, 
ciliates and most excavates, a slightly more restricted set 
of proteins was observed (Table 1). Notably, none of 
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them harboured the Apc5 subunit of the structural 
complex, along with two components of the TPR arm 
(Apcl2 and Apcl3), whereas we detected Apc4 only in 
the ciliate Tetrahymena thermophila and the haptophyte 
E. huxleyi, and Apc7 in the excavate Trichomonas vagi- 
nalis (Table 1). This suggested that the great similarity 
observed in the distribution of APC/C components in 
these four major eukaryotic lineages was likely due to 
convergent losses during their evolution (Figure 2). The 
situation was completely different in Apicomplexa, the 
second main lineage of alveolates in our dataset together 
with ciliates. While the three ciliates have retained a 
large number of components inferred to be present in 
LECA (Table 1), most of them have been lost in Api- 
complexa (Table 1). More specifically, all APC/C pro- 
teins were missing in Babesia bovis and Theileria 
annulata, whereas the remaining four apicomplexan 
species harboured only four or five of them. Surpris- 
ingly, components present in those organisms were 
diverse depending on the species (Table 1). For instance, 
ApclO, Apcll and Apcl were found in Toxoplasma 
gondii, whereas the two Plasmodium species contained 
orthologues of ApclO, Apcll, Apc3 and the anaerobic 
Cryptosporidium hominis harboured Apcl, Apcll, Apc3, 
Apc6 and Apc8 (Additional file 1, Table SI). The pre- 
sence of nearly all components in ciliates indicated that 
massive and differential losses occurred secondarily in 
Apicomplexa (Figure 2). As mentioned above, such mas- 
sive losses were also observed in the excavate G. intesti- 
nalis and in the amoebozoan E. histolytica (Additional 
file 1, Table SI). However, it is important to note that 
when orthologues existed in these lineages, they showed 
highly divergent sequences compared to those found in 
other eukaryotic lineages (not shown). It is thus possible 
that orthologues of some components might have 
escaped detection (despite the use of PSI-BLAST 
searches) because of their extreme degree of sequence 
divergence. In any case, the possible massive losses or 
the high divergence of APC/C components both sug- 
gested that important changes have occurred relatively 
recently in these parasitic lineages, maybe as a conse- 
quence of their atypical cell division mechanism. This is 
notably the case of Theileria that, acting like a disguised 
chromosome during host cell division, inserts itself into 
both daughter cells by co-opting parts of the host cell 
division machinery, in particular the host cell's microtu- 
bules to be pulled towards the opposing ends of the 
dividing cell [25]. 

An interesting evolutionary pattern emerged from our 
analyses concerning the APC/C adaptors/co-activators. 
Their phylogenies supported that only the two paralo- 
gues Cdc20 and Cdhl were present in LECA and con- 
served in nearly all eukaryotic lineages (even in 
Microsporidia), whereas all the remaining co-activators 



resulted from independent duplications that occurred 
recently in different eukaryotic lineages (Additional file 

4, Figure SI and Additional file 5, Figure S2). For 
instance, Rap and Cortex resulted from duplications of 
Cdhl and Cdc20, respectively, which occurred in D. 
melanogaster, whereas Amal and Mfrl derived from 
duplications of Cdc20 in Ascomycota and Cdhl in S. 
pombe, respectively. Within plants, our analyses con- 
firmed the presence of multiple Cdc20 and Cdhl copies 
in A. thaliana, but also in other land plants (Oryza 
sativa, Selaginella moellendorfii and Physcomitrella 
patens) (Additional file 1, Table SI) [19]. Phylogenetic 
analyses suggested that the numerous homologues 
observed in these plants arose from independent dupli- 
cations (Additional file 4, Figure SI and Additional file 

5, Figure S2). In A. thaliana, the various Cdc20 and 
Cdhl paralogues have been shown to be differently 
expressed through the cell cycle and depending on cell 
types or tissues [19], suggesting that subfunctionaliza- 
tion events occurred after the duplications. More intri- 
guing was the huge expansion of the repertory of 
adaptor/co-activators observed in the two ciliates T. 
thermophila and Paramecium tetraurelia, for which we 
identified eight and ten copies of Cdhl, respectively. 
Among excavates, Leishmania and Trypanosoma gen- 
omes encoded only one adaptor/co-activator affiliated to 
the Cdc20 subfamily (Table 1), whereas the genome of 
Naegleria (as in most eukaryotes) encoded one Cdc20 
and one Cdhl copies. The genome of Trichomonas con- 
tained three homologues but due to their great diver- 
gence we were unable to classify them as Cdc20 or 
Cdhl without ambiguity (Additional file 1, Table SI). 
Finally, in G. intestinalis as in some apicomplexa (Babe- 
sia, Theileria and Toxoplasma), we failed to detect any 
adaptor/co-activator, reinforcing the hypothesis that 
their APC/C proteins have experienced a very divergent 
and fast evolution. 

Regarding the main APC/C targets, our phylogenetic 
analyses were not conclusive in the case of cyclins A 
and B and Cdks-1 and -2 to determine whether they 
were found in LECA or not because these proteins 
belonged to very large multigenic families with complex 
evolutionary histories precluding a precise inference of 
their evolutionary origin (not shown). For the remaining 
targets, our analyses allowed inferring that the separase 
and the nine subunits composing the CC (i.e. the core 
complex, associated proteins, the loading complex and 
the proteins involved in the cohesion establishment) 
were present in LECA (Table 2 and Additional file 2, 
Table S2) and have been conserved in most eukaryotic 
lineages (Table 2 and Figure 2). The taxonomic distribu- 
tion of Smcl, Smc3, Sccl, Scc3 and Psd5 homologues 
was globally in agreement with a previous study focused 
on the analysis of 29 genes involved in meiosis in 
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eukaryotes [26]. In contrast, Scc4 (a protein composing 
the loading complex) and Wpll/Rad61 (a protein that 
forms together with Psd5 a complex involved in CC 
maintenance [12,13]) were present only in Viridiplan- 
tae and in Opisthokonta suggesting convergent losses 
in other lineages (Table 2 and Figure 2). However, it 
can also be speculated that these proteins are not 
under strong selective pressure, as attested by the fast 
evolutionary rate of Wpll/Rad61 found in Saccharo- 
mycetaceae that are shorter and highly divergent com- 
pared to those of metazoa (647 versus 1190 amino 
acids, 14% identity and 29% similarity between yeast 
(NP_010297) and human (NP_055860) sequences). So 
it would even be possible that they have been replaced 
by non-homologous proteins in other lineages. The 
only targets of APC/C that were not inferred to be 
present in LECA are securins that were found only in 
Metazoa and Fungi. However, even though they fulfil 
the same function (preventing the separation of the 
two sister chromatids) through the binding of 
separases that are homologous in metazoan and fungal 
species, fungal securins are not homologous to those 
from metazoa [27], suggesting again a non-homologous 
replacement in one of these two groups (Table 1 and 
Figure 2). 

Functional data point to a nearly modern APC/C 
controlling the cell cycle in LECA 

Our phylogenomic analysis of the APC/C, its main 
adaptors/co-activators and targets supported the hypoth- 
esis that most of the corresponding genes were already 
present in LECA. It was therefore tempting to conclude 
that a complex similar to the one involved in present- 
day regulation of the eukaryotic cell cycle was already 
present in this ancestral organism. Supporting this 
hypothesis, nearly all the proteins under study possessed 
conserved functional domains similar to those present 
in components that were experimentally proved to be 
part of the APC/C or targeted by this ubiquitin ligase 
(Additional file 6, Table S3 and Additional file 7, Table 
S4). This indicated that all the identified orthologues of 
APC/C components and targets had similar molecular 
functions. Moreover, experimental data from plants and 
excavates have shown that most components identified 
by our analysis were part of or targeted by the APC/C 
[18-21,28-35]. This strongly suggested that they inher- 
ited their function from the ancestral proteins present 
in LECA, and therefore that a nearly modern APC/C 
and control of the cell cycle existed in LECA. This also 
suggested that, in lineages for which no functional data 
were available, the orthologues that we have identified 
were likely involved in the control of the cell cycle and, 
therefore, may constitute interesting targets for experi- 
mental work. 



Although the origin of most APC/C components and 
targets could be traced back to LECA and most of them 
have been conserved throughout the evolution of eukar- 
yotes, some component gains could also be observed. 
Most of them resulted from gene duplications of adap- 
tors/co-activators that occurred independently in differ- 
ent lineages. This was in agreement with recent reports 
of new, and often specific, activities of the APC/C in 
some eukaryotes (e.g. the neuronal activity [36] or the 
early eye development [16], or the progression of the 
endocycle in flowering plant endosperm [18,20]). This 
suggested that most of the APC/C evolution since 
LECA has concerned the acquisition of new regulatory 
functions by increasing the repertoire of adaptors/co- 
activators, even if we can not rule out the possibility 
that adaptors/co-activators present in single copies in 
some lineages were (and still are) able to interact with a 
larger spectrum of targets than their multiple-copies 
counterparts. In addition to the classical activators/co- 
activators, a recent interactomic study in A. thaliana 
suggested the presence of three novel APC/C interactors 
specific to land plants that were not homologous to the 
Cdc20/Cdhl family [34]. However, although they inter- 
acted with the APC/C, their biological function (e.g. 
adaptor/co-activators, inhibitors, targets, etc.) has still 
not been established. Nevertheless, this supported that 
lineage-specific innovations are expected to be discov- 
ered when biological data on a broader diversity of 
eukaryotes becomes available. 

In contrast, we also observed that convergent events 
of streamlining occurred secondarily in various lineages, 
like Apicomplexa, G. lamblia and E. histolytica. The 
reasons explaining those massive loss events are unclear, 
though we could not discard that, at least for a number 
of cases, they might be linked to an extreme accelera- 
tion of their evolutionary rate beyond detection. This 
hypothesis was supported by the high evolutionary rates 
exhibited by the few components still harboured by Api- 
complexa. However, the detection of APC/C subunits 
and targets in fast evolving organisms, like E. cuniculi 
and T. vaginalis, rather suggested that most of the miss- 
ing components in G. lamblia, E. histolytica and Api- 
complexa, reflected true losses. In that case we could 
wonder whether the lost components have been 
replaced by non homologous proteins that fulfil the 
same role or whether these parasites are able to recruit 
the APC/C components from their hosts. In both cases, 
experimental investigations in these parasitic lineages 
will be useful to elucidate the nature of their APC/C or 
even to discover putative divergent systems involved in 
the control of the cell cycle that may provide interesting 
medical drug targets. Likewise, a previous phylogenomic 
study of proteins involved in late cytokinesis revealed a 
similar pattern of reductive evolution in these lineages 
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[37]. Indeed, E. histolytica, G. lamblia and Apicomplexa 
have undergone massive losses of proteins of the cyto- 
kinesis machinery, including conserved ancient ones 
inferred to have been present in LECA [37]. This infor- 
mation combined to our present analysis suggests that 
major changes have occurred in various steps of the cell 
cycle in these parasitic eukaryotes. 

Most components of the APC/C and targets are 
eukaryotic innovations 

Despite our extensive survey of public sequence data- 
bases, we did not identify any homologue of the APC/C 
components in prokaryotes with two exceptions. This 
indicated that this large E3 complex and its main targets 
are eukaryotic innovations that emerged after the 
separation of this domain from prokaryotes but prior to 
its diversification into the present-day eukaryotic 
lineages. The two exceptions, Smcl and Smc3, are two 
paralogous proteins that are part of the core complex of 
the cohesin complex (Table 1). According to previous 
reports and to their critical role in higher-order chro- 
mosome organization and dynamics [38,39], we identi- 
fied homologues of Smcl and Smc3 in nearly all 
archaeal and bacterial lineages (not shown). The lack of 
APC/C prokaryotic homologues was surprising because 
distant homologues harbouring structures similar to 
eukaryotic ubiquitin, El and E2 exist in prokaryotes [40] 
and because a bona fide homologue of the eukaryotic 
proteasome has been described in Archaea and Actino- 
bacteria [41,42]. Moreover, it was recently reported that 
homologues of the eukaryotic ubiquitination pathway 
are encoded in the genome of the archaeon 'Candidatus 
Caldiarchaeum subterraneum' [43], a relative of the 
recently proposed phylum Thaumarchaeota [44]. This 
system is composed of a cluster of four genes coding for 
the ubiquitin, El -like and E2-like enzymes and a small 
Zn RING finger protein [43]. The first three proteins 
are much more similar to their eukaryotic counterparts 
than to the very distant homologues usually found in 
prokaryotes. Since no bona fide homologue of E3 
enzymes has been identified in this archaeon, it was pro- 
posed that the fourth protein might mediate the ligation 
of ubiquitin [43]. Despite this recent discovery, our ana- 
lyses together with data from the literature suggested 
that the vast majority of APC/C components and main 
targets were eukaryotic innovations that may have 
played a role in the emergence and evolution of the 
complex cell cycle observed in this domain. 

APC/C components and main targets can be used to infer 
the phylogeny of eukaryotes 

Proteins inferred to have been present in LECA are 
valuable material to reconstruct the characteristics of 
this ancestral organism (e.g. [37,45-49]). In addition, 



they can preserve a phylogenetic signal useful to infer 
the evolutionary history of eukaryotes. Until now, most 
analyses dedicated to the reconstruction of the eukaryo- 
tic phylogeny were based on the analysis of components 
of informational systems (those involved in the trans- 
mission and expression of the genetic information). This 
was so because most of the genes coding for these pro- 
teins present the advantage of being (i) part of mono- 
genic gene families allowing the easy identification of 
orthologous proteins, (ii) slowly evolving, (iii) ancient 
and well conserved among life domains, and (iv) rarely 
exchanged by horizontal gene transfers. Accordingly, 
they represent first choice material to investigate ancient 
evolution in all domains of life [50-53]. Phylogenetic 
studies of the eukaryotic domain did not escape this 
rule and most of them have been largely based on the 
phylogenetic analysis of informational proteins (e.g., 
[51,54-57]). By contrast, most proteins involved in 
housekeeping functions (operational proteins) are con- 
sidered to evolve faster than those of informational sys- 
tems, and thus to be less suitable to study ancient 
evolution. Moreover, they are often part of large and 
complex protein families that have experienced numer- 
ous gene duplication and loss events during their evolu- 
tionary history, meaning that distinguishing between 
orthologues and paralogues is difficult and requires fas- 
tidious preliminary analyses. Consequently, although 
these proteins are more numerous than informational 
ones and often of larger size, they have rarely been used 
to infer the ancient evolution of eukaryotes. Neverthe- 
less, typical datasets based on informational proteins 
have been shown to be insufficient to robustly infer all 
the deep nodes of the global eukaryotic phylogenies 
[58,59]. Increasing the protein sampling is therefore 
becoming as necessary as increasing the taxonomic sam- 
pling in order to fully resolve the phylogeny of eukar- 
yotes, meaning that new useful protein markers have to 
be found among the conserved operational proteins. 

Our phylogenetic analyses of the APC/C subunits and 
main targets showed that with a few exceptions they are 
ancient proteins well conserved throughout the diversity 
of present day eukaryotic lineages. Accordingly, they are 
potential suitable markers to reconstruct global eukaryo- 
tic phylogenies. The maximum likelihood and Bayesian 
phylogenetic analyses of individual components showed 
that they have retained ancient phylogenetic signal 
despite the fact that some basal nodes of the inferred 
phylogenies showed a poor resolution (a typical result in 
single-marker analyses) (Additional file 3, Data SI). To 
increase resolution, we built a supermatrix (3115 posi- 
tions) by concatenating the individual alignments of the 
24 APC/C components and main targets inferred to be 
present in LECA. To minimise potential long-branch 
attraction (LBA) tree reconstruction artefacts [60], we 
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discarded at this step the extremely fast evolving and/or 
very partial sequences of Enterocytozoon bieneusi 
(Microsporidia), Plasmodium spp. (Apicomplexa) and B. 
hominis (Heterokonta), all them having > 60% of miss- 
ing data. We also discarded T. vaginalis (Excavata) 
because this species contained multiple (and often diver- 
gent) copies of several APC/C components and main 
targets. The resulting Bayesian and maximum likelihood 
trees recovered the monophyly of each eukaryotic phyla 
and most often supported by high Bayesian posterior 
probabilities (PP) and maximum likelihood bootstrap 
values (BV): Discicristata (PP = 0.98; BV = 97%), Het- 
erokonta (PP = 1.0; BV = 100%), Metazoa (PP = 1.0; BV 
= 100%), Fungi (PP = 0.87; BV = 98%), Choanoflagellata 
(PP = 1.0; BV = 100%), Viridiplantae (PP = 1.0; BV = 
99%) and Alveolata (PP = 0.95 and BV = 78%) (Figure 
3). The most remarkable result was the good support 
retrieved for the monophyly of Fungi, including the 
extremely fast evolving Microsporidia. The APC/C pro- 
teins are therefore among the rare ones [61] that sup- 
port the correct phylogenetic placement of these 
parasites as members of the Fungi. 

Regarding the relationships among the main eukaryo- 
tic groups, the sister-grouping of the two choanoflagel- 
lates and C. owczarzaki was supported both by the 
Bayesian (PP = 0.88) and maximum likelihood (BV = 
87%) trees; these two groups being closely related to 
Metazoa (namely, a support for the Holozoa with PP = 
0.97 and BV = 100%). This was in agreement with pre- 
vious results [62,63], but not with those supporting the 
emergence of Capsaspora at the base of the group com- 
posed of Choanoflagellata and Metazoa [64,65]. The 
phylogenetic signal contained in APC/C components 
and main targets also supported the monophyly of 
Opisthokonta (PP = 0.91; BV = 100%). Moreover, if we 
accept the putative unikont-bikont rooting [23], our 
analyses also retrieved the Unikonta, including the apu- 
sozoan Thecamonas trahens, albeit with a weaker sup- 
port (PP = 0.81; BV = 78%). Similarly, most of the 
relationships among bikont lineages were poorly 
resolved (Figure 3). In particular, the monophyly of 
Plantae was not recovered: green plants and algae 
appeared closely related to the haptophyte E. huxleyi 
(PP = 0.78; BV = 81%), whereas the phylogenetic posi- 
tion of red algae was not resolved. This agreed with 
other studies showing the difficulty to retrieve the 
monophyly of Plantae, even with much larger datasets 
(e.g., [66]). Finally, heterokonts formed the sister group 
of green plants and haptophytes (PP = 0.89; BV = 
100%). 

Taken together, our phylogenetic analyses showed that 
despite the relative small size of the supermatrix con- 
structed with the APC/C components and main targets 
(3115 positions), it contained an interesting phylogenetic 



signal that can be useful to infer part of the ancient evo- 
lutionary history of eukaryotes. However, the large size 
of most of the APC/C subunits and main targets made 
difficult the use of the abundant body of sequence data 
coming from the analysis of expressed sequence tags 
(ESTs) because of their incomplete sequence coverage, 
which limited our taxonomic sampling. Nevertheless, 
our dataset can now be employed as a reference to 
enrich it with assembled EST sequences as well as with 
those from upcoming complete genome sequences from 
various protist species. This limitation will likely apply 
to other operational proteins candidates to be used as 
markers for the reconstruction of global eukaryotic phy- 
logenies but, as already discussed, their inclusion in 
multi-marker analyses becomes crucial to infer several 
deep nodes that the traditional supermatrices of infor- 
mational proteins have failed to resolve. 

Conclusions 

Our analyses of the APC/C and its main targets showed 
that this complex system was very likely present in 
LECA and has been conserved, to a few exceptions, all 
along the diversification of the eukaryotic domain. This 
study provided first insights into the mechanisms 
responsible of the control of the cell cycle in LECA, sug- 
gesting that it was tightly regulated like in present-day 
eukaryotes. Finally we showed that the components of 
the APC/C and its main targets can be good phyloge- 
netic markers to complement those used so far (mostly 
proteins involved in transcription and translation). 
Indeed, the latter have proven not to be sufficient to 
fully resolve the phylogeny of eukaryotes, making neces- 
sary to identify new complementary markers. This will 
certainly be a difficult task that will require many ana- 
lyses but we think that the phylogenomic study of con- 
served cellular systems is a promising approach to 
tackle this issue. 

Methods 

Dataset assembly 

We used the 37 APC/C components and main targets 
identified in four opisthokont species {Homo sapiens, D. 
melanogaster, S. cerevisiae and 5. pombe [7,11]), the 
plant A. thaliana [17,19] and the kinetoplastid T. brucei 
[21] (Table 1) to survey public sequence databases. We 
identified homologues of these proteins using BLASTp 
and PSI-BLAST (with default parameters) [67] in a sub- 
set of complete or ongoing genomes representative of 
eukaryotic diversity available at the NCBI (Figure 2). To 
increase the taxonomic sampling, homologues of Mono- 
siga brevicollis, Salpingoeca rosetta, Lottia gigantea, 
Nematostella vectensis, Helobdella robusta, Daphnia 
pulex, Capsaspora owczarzaki, Batrachochytrium den- 
drobatidis, Spizellomyces punctatus, Thecamonas 



Erne et al. BMC Evolutionary Biology 201 1 , 11 :265 
http://www.biomedcentral.eom/1 471 -2 1 48/1 1 /265 



Page 11 of 14 



■ Nosema ceranae ~\ 



.91 



97- 



.81 



IT 



■ Spizellomyces punctatus 

Batrachochytrium dendrobatidis 

■ Ustilago maydis 



£■ Encephalitozoon intestinalis 
"7*— Encephalitozoon cuniculi 



A 



■ Cryptococcus neoformans 

■ Schizosaccharomyces pombe 
■ Saccharomyces cerevisiae 

■ Neurospora crassa 



Aspergillus fumigatus 
Capsaspora owczarzaki — 
— Salpingoeca rosettt 

- Monosiga brevicollis . 



.95 

.9/ 



Trichoplax adhaerens 

Nematostella vectensis 

r — Branchiostoma floridae 
IT c Homo sapiens 
<t- r 



.97'' 



.95 
.76 



f" - Danio rerio 
e — Lottia gigantea 
E Helobdella robusta 



■ Caenorhabditis elegans 



■ Brugia malayi 



7«l 



— Drosophila melanogaster 
■ Daphnia pulex 



78 1 — Apis mellifera 



■ Dictyostelium discoideum 



Thecamonas trahens 



■ Naegleria gruberi 




Xl* — 7 
^ # 



Trypanosoma cruzi 
Trypanosoma brucei 
Leishmania major 
Leishmania infantum 



Galdieria sulphuraria ' 



■ Cyanidioschyzon merolae . 



Oxytricha trifallax 

£ Tetrahymena thermophila 

i* Paramecium tetraurelia 



■ Cryptosporidium hominis 



Phytophthora ramorum 
>* Phytophthora infestans 

Thalassiosira pseudonana 
Phaeodactylum tricornutum 
Aureococcus anophagefferens 

Ectocarpus siliculosus 

Emiliania huxleyi 



Selaginella moellendorffii 

Physcomitrella patens 
Oryza sativa 
Arabidopsis thaliana 
Micromonas pusilla 
Ostreococcus tauri 
99^- Ostreococcus lucimarinus 
Chlorella vulgaris 
Volvox carter! 

- Chlamydomonas reinhardtii . 



Fungi 



Capsaspora 



^ Choanoflagellata 



Metazoa 



- Amoebozoa 

- Apusozoa 

Discicristata 



J) Rhodophyta 



Alveolata 



Heterokonta 



Haptophyta 



Viridiplantae 



Figure 3 Bayesian phylogenetic tree of eukaryotes. The tree was inferred by PhyloBayes based on a supermatrix made of the 24 APC/C 
components, adaptors/co-activators and main targets inferred to be present in the last common ancestor of eukaryotes (3115 amino-acid 
positions, 56 species). Numbers at nodes are Bayesian posterior probabilities, filled circles indicate maximum likelihood bootstrap values of 100% 
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trahens, Naegleria gruberi, Phaeodactylum tricornutum, 
Aureococcus anophagefferens, Ostreococcus lucimarinus, 
Physcomitrella patens, Chlorella vulgaris, Micromonas 
pusilla, Selaginella moellendorffii and Emiliania huxleyi 
were retrieved using the BLASTp and tBLASTn pro- 
grams from the JGI http://genome.jgi-psf.org/, the Broad 
Institute http://www.broadinstitute.org and the TBestDB 
database http://tbestdb.bcm.umontreal.ca/ [68]. In addi- 
tion, homologues of two representatives of Rhodophyta 



were retrieved from the Galdieria sulphuraria genome 
project http://genomics.msu.edu/cgi-bin/galdieria/blast. 
cgi and the Cyanidioschyzon merolae genome project 
http://merolae.bioLs.u-tokyo.ac.jp/blast/blast.html. 
BLAST outputs were examined by eye to identify homo- 
logues of each protein to avoid applying an arbitrary 
cutoff on e-value or score. To ensure an exhaustive 
sampling of homologues, we performed additional 
searches using as seeds homologues that were identified 
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at previous steps. The absence of any homologue in a 
given lineage was systematically verified by hand using 
tBLASTn searches on the nucleotide sequences of the 
corresponding complete genomes. 

For each protein, the homologous sequences were 
gathered in a dataset and aligned with MAFFT 6.833 
[69]. All the resulting alignments were edited and manu- 
ally adjusted using the ED program from the MUST 
package [70]. Regions where homology between sites 
was doubtful were manually removed from the align- 
ments before phylogenetic analyses. A supermatrix was 
built by concatenating the alignments corresponding to 
the 24 APC/C components, adaptor/co-activators and 
targets inferred to be present in LECA. 

Conserved functional domain search 

The identification of functional domains was carried out 
using the HMMER package [71] and the HMM profiles 
of the Pfam database [72]. HMM profiles having e- 
values lower than 0.1 were considered as significant. 

Phylogenetic analysis 

Phylogenetic trees were reconstructed for each single 
protein dataset, except for Apc9, Apcl4, Apcl5, Apcl6 
and securin because of their very restricted number of 
homologues. Maximum Likelihood (ML) phylogenetic 
trees were computed using Treefinder [73] with the LG 
model [74] and a gamma correction to take into account 
the heterogeneity of evolutionary rates across sites (4 
discrete classes of sites and estimated alpha parameter), 
as proposed by the model selection tool implemented in 
Treefinder. The branch robustness of ML trees was esti- 
mated using the non parametric bootstrap procedure 
(100 replicates of the original alignment) implemented 
in Treefinder with the same parameters than for ML 
tree inference. Bayesian trees were computed using 
MrBAYES 3.1.2 [75] with a mixed amino acid substitu- 
tion model and a gamma correction (four rate cate- 
gories). The Markov chain Monte Carlo search was run 
with 4 chains for 1, 000, 000 generations, with trees 
being sampled every 100 generations (the first 2, 500 
trees were discarded as "burn-in"). Individual Bayesian 
trees are provided as Additional file 3, Data SI. 

TreeFinder and PhyloBayes 3.2 [76] were used to per- 
form ML and Bayesian analyses on the concatenated 
dataset. Treefinder was run with the same parameters 
than for the analysis of single datasets. Phylobayes was 
run using the CAT model and a gamma correction with 
4 rate categories [77]. For each dataset, two independent 
chains were run for at least 10000 cycles, saving one 
tree in ten. The first 300 trees were discarded as "burn- 
in", and the remaining trees from each chain in each 
dataset were used to test for convergence and compute 
the 50% majority rule consensus tree. 



Inference of the origin of APC/C components and main 
targets 

The phylogenetic analysis of each individual component 
allowed distinguishing orthologues that originated from 
speciation events from paralogues that resulted from 
gene duplications. This point is important because infer- 
ring the evolutionary origin of a protein requires the 
analysis of the evolutionary history of orthologues. To 
determine the origin of each component and subunit, 
we used a parsimony criterion, making the assumption 
that horizontal gene transfers between eukaryotes are 
rare. Accordingly, the ancestral absence of a component 
in the ancestor of a taxonomic group was inferred if no 
orthologues are found in any present-day representative 
of the group (even in those that were not included in 
our taxonomic sampling). By contrast, the presence of 
orthologues in representatives of two lineages was inter- 
preted as the existence of the corresponding component 
in their last common ancestor. Similarly, the presence of 
a component in an ancestral taxon was also inferred if 
orthologues are present in at least one representative of 
its offspring and in one of its ancestors. 

Additional material 



Additional file 1: Table SI. Table showing the taxonomic distribution of 
homologues of APC/C subunits and activators in eukaryotes. 

Additional file 2: Table S2. Table showing the taxonomic distribution of 
homologues of APC/C main targets in eukaryotes. 

Additional file 3: Data SI. Bayesian phylogenetic trees inferred MrBayes 
(*.con files) and ML phylogenetic trees inferred with TreeFinder (*.tre files 
= ML topologies; *.boot.trees = topologies generated by the bootstrap 
analysis) for each APC/C subunit, activator and main target (except Apc9, 
Apc14, Apc15 and Apc16 due to a too restricted number of 
homologues). 

Additional file 4: Figure SI. Bayesian phylogenetic tree of the adaptor/ 
co-activator Cdhl inferred with Mrbayes. Numbers at nodes correspond 
to posterior probabilities (PP) estimated by MrBayes and bootstrap values 
(BV) estimated by TreeFinder (only PP and BV greater than 0.5 and 50%, 
respectively, are shown). The scale bar indicates the average number of 
substitutions per site. 

Additional file 5: Figure S2. Bayesian phylogenetic tree of the adaptor/ 
co-activator Cdc20 inferred with Mrbayes. Numbers at nodes correspond 
to posterior probabilities (PP) estimated by MrBayes and bootstrap values 
(BV) estimated by TreeFinder (only PP and BV greater than 0.5 and 50%, 
respectively, are shown). The scale bar indicates the average number of 
substitutions per site. 

Additional file 6: Table S3. Table showing the conserved functional 
domains present in homologues of APC/C subunits and activators. 

Additional file 7: Table S4. Table showing the conserved functional 
domains present in homologues of APC/C main targets. 



Acknowledgements and Funding 

CBA is supported by an Action Thematique et Incitative sur Programme 
(ATIP) of the French Centre National de la Recherche Scientifique (CNRS). We 
would like to thank the JGI institute, the TBestDB database, the Galdieria 
sulphuraria genome project and the Cyanidioschyzon merolae genome 
project for sharing sequence data from ongoing projects with the scientific 



Erne ef al. BMC Evolutionary Biology 201 1 , 11 :265 
http://www.biomedcentral.eom/1 471 -2 1 48/1 1 /265 



Page 13 of 14 



community. We also thank Rym Agrebi and Celine Petitjean for careful 
reading of the manuscript and stimulating discussions. 

Author details 

1 Aix-Marseille Universite, Laboratoire de Chimie Bacterienne - UPR CNRS 
9043, Marseille, France. 2 Unite d'Ecologie, Systematique et Evolution - UMR 
CNRS 8079, Universite Paris-Sud, Orsay, France. 3 Present Address: Universite 
de Lyon; Universite Lyon 1; CNRS; UMR5558, Laboratoire de Biometrie et 
Biologie Evolutive, 43 boulevard du 1 1 novembre 1918, Villeurbanne, France. 

Authors' contributions 

LE, DM and CBA designed the analyses and wrote the paper. LE and AT 
performed the analyses. All authors read and approved the final manuscript. 

Received: 18 April 2011 Accepted: 23 September 2011 
Published: 23 September 2011 

References 

1 . Thornton BR, Toczyski DP: Precise destruction: an emerging picture of the 
APC. Genes Dev 2006, 20(22)3069-3078. 

2. Fang S, Weissman AM: A field guide to ubiquitylation. Cell Mol Life Sci 
2004, 61 (13):1 546-1 561. 

3. Pickart CM: Mechanisms underlying ubiquitination. Annu Rev Biochem 
2001, 70:503-533. 

4. Schreiber A, Stengel F, Zhang Z, Enchev Rl, Kong EH, Morris EP, 
Robinson CV, da Fonseca PC, Barford D: Structural basis for the subunit 
assembly of the anaphase-promoting complex. Nature 201 1, 
470(7333)227-232. 

5. Reed SI: Ratchets and clocks: the cell cycle, ubiquitylation and protein 
turnover. Nor Rev Mol Cell Biol 2003, 4(1 1):855-864. 

6. Manchado E, Eguren M, Malumbres M: The anaphase-promoting complex/ 
cyclosome (APC/C): cell-cycle-dependent and -independent functions. 
Biochem Soc Trans 2010, 38(Pt 1):65-71. 

7. Peters JM: The anaphase promoting complex/cyclosome: a machine 
designed to destroy. Nat Rev Mol Cell Biol 2006, 7(9):644-656. 

8. Yoon HJ, Feoktistova A, Wolfe BA, Jennings JL, Link AJ, Gould KL: 
Proteomics analysis identifies new components of the fission and 
budding yeast anaphase-promoting complexes. Curr Biol 2002, 
12(23)2048-2054. 

9. Hutchins JR, Toyoda Y, Hegemann B, Poser I, Heriche JK, Sykora MM, 
Augsburg M, Hudecz O, Buschhorn BA, Bulkescher J, ef al: Systematic 
analysis of human protein complexes identifies chromosome 
segregation proteins. Science 2010, 328(5978)593-599. 

10. Acquaviva C, Pines J: The anaphase-promoting complex/cyclosome: APC/ 
C. J Cell Sci 2006, 1 19(Pt 12)2401-2404. 

1 1. Matyskiela ME, Rodrigo-Brenni MC, Morgan DO: Mechanisms of ubiquitin 
transfer by the anaphase-promoting complex. J Biol 2009, 8(10):92. 

12. Nasmyth K, Haering CH: Cohesin: its roles and mechanisms. Annu Rev 
Genet 2009, 43:525-558. 

13. Onn I, Heidinger-Pauli JM, Guacci V, Unal E, Koshland DE: Sister chromatid 
cohesion: a simple concept with a complex reality. Annu Rev Cell Dev Biol 
2008, 24:105-129. 

14. Oelschlaegel T, Schwickart M, Matos J, Bogdanova A, Camasses A, Havlis J, 
Shevchenko A, Zachariae W: The yeast APC/C subunit Mnd2 prevents 
premature sister chromatid separation triggered by the meiosis-specific 
APC/C-Amal. Cell 2005, 120(6):773-788. 

1 5. Pesin JA, Orr-Weaver TL: Developmental role and regulation of cortex, a 
meiosis-specific anaphase-promoting complex/cyclosome activator. PLoS 
Genet 2007, 3(1 1):e202. 

16. Pimentel AC, Venkatesh TR: rap gene encodes Fizzy-related protein (Fzr) 
and regulates cell proliferation and pattern formation in the developing 
Drosophila eye-antennal disc. Dev Biol 2005, 285(2):436-446. 

1 7. Capron A, Okresz L, Genschik P: First glance at the plant APC/C, a highly 
conserved ubiquitin-protein ligase. Trends Plant Sci 2003, 8(2):83-89. 

18. Eloy NB, Coppens F, Beemster GT, Hemerly AS, Ferreira PC: The Arabidopsis 
anaphase promoting complex (APC): regulation through subunit 
availability in plant tissues. Cell Cycle 2006, 5(1 7):1 957-1 965. 

19. Lima Mde F, Eloy NB, Pegoraro C, Sagit R, Rojas C, Bretz T, Vargas L, 
Elofsson A, de Oliveira AC, Hemerly AS, ef al: Genomic evolution and 
complexity of the Anaphase-promoting Complex (APC) in land plants. 
BMC Plant Biol 2010, 10:254. 



20. Fulop K, Tarayre S, Kelemen Z, Horvath G, Kevei Z, Nikovics K, Bako L, 
Brown S, Kondorosi A, Kondorosi E: Arabidopsis anaphase-promoting 
complexes: multiple activators and wide range of substrates might keep 
APC perpetually busy. Cell Cycle 2005, 4(8):1 084-1 092. 

21. Kumar P, Wang CC: Depletion of anaphase-promoting complex or 
cyclosome (APC/C) subunit homolog APC1 or CDC27 of Trypanosoma 
brucei arrests the procyclic form in metaphase but the bloodstream 
form in anaphase. J Biol Chem 2005, 280(36)31783-31791. 

22. Eisen JA, Fraser CM: Phylogenomics: intersection of evolution and 
genomics. Science 2003, 300(5626)4 706-1707. 

23. Stechmann A, Cavalier-Smith T: The root of the eukaryote tree 
pinpointed. Curr Biol 2003, 1 3(1 7):R665-666. 

24. Seidl MF, Schultz J: Evolutionary flexibility of protein complexes. BMC Evol 
Biol 2009, 9:155. 

25. von Schubert C, Xue G, Schmuckli-Maurer J, Woods KL, Nigg EA, 
Dobbelaere DA: The transforming parasite Theileria co-opts host cell 
mitotic and central spindles to persist in continuously dividing cells. 
PLoS Biol 2010, 8(9). 

26. Malik SB, Pightling AW, Stefaniak LM, Schurko AM, Logsdon JM Jr: An 
expanded inventory of conserved meiotic genes provides evidence for 
sex in Trichomonas vaginalis. PLoS One 2008, 3(8):e2879. 

27. Zou H, McGarry TJ, Bernal T, Kirschner MW: Identification of a vertebrate 
sister-chromatid separation inhibitor involved in transformation and 
tumorigenesis. Science 1999, 285(5426):41 8-422. 

28. Lam WS, Yang X, Makaroff CA: Characterization of Arabidopsis thaliana 
SMC1 and SMC3: evidence that AtSMC3 may function beyond 
chromosome cohesion. J Cell Sci 2005, 1 18(Pt 14)3037-3048. 

29. Li W, Yang X, Lin Z, Timofejeva L, Xiao R, Makaroff CA, Ma H: The 
AtRAD51C gene is required for normal meiotic chromosome synapsis 
and double-stranded break repair in Arabidopsis. Plant Physiol 2005, 
138(2)365-976. 

30. Liu Cm CM, McElver J, Tzafrir I, Joosen R, Wittich P, Patton D, Van 
Lammeren AA, Meinke D: Condensin and cohesin knockouts in 
Arabidopsis exhibit a titan seed phenotype. Plant J 2002, 29(4):405-415. 

31. Bai X, Peirson BN, Dong F, Xue C, Makaroff CA: Isolation and 
characterization of SYN1, a RAD21 -like gene essential for meiosis in 
Arabidopsis. Plant Cell 1999, 11(3):41 7-430. 

32. Bhatt AM, Lister C, Page T, Fransz P, Findlay K, Jones GH, Dickinson HG, 
Dean C: The DIF1 gene of Arabidopsis is required for meiotic 
chromosome segregation and belongs to the REC8/RAD21 cohesin gene 
family. Plant J 1999, 1 9(4):463-472. 

33. Cai X, Dong F, Edelmann RE, Makaroff CA: The Arabidopsis SYN1 cohesin 
protein is required for sister chromatid arm cohesion and homologous 
chromosome pairing. J Cell Sci 2003, 1 16(Pt 14)2999-3007. 

34. Van Leene J, Hollunder J, Eeckhout D, Persiau G, Van De Slijke E, Stals H, 
Van Isterdael G, Verkest A, Neirynck S, Buffel Y, et al: Targeted 
interactomics reveals a complex core cell cycle machinery in 
Arabidopsis thaliana. Mol Syst Biol 2010, 6397. 

35. Sebastian J, Ravi M, Andreuzza S, Panoli AP, Marimuthu MP, Siddiqi I: The 
plant adherin AtSCC2 is required for embryogenesis and sister- 
chromatid cohesion during meiosis in Arabidopsis. Plant J 2009, 
59(1):1-13. 

36. Harper JW, Burton JL, Solomon MJ: The anaphase-promoting complex: it's 
not just for mitosis any more. Genes Dev 2002, 16(17)2179-2206. 

37. Erne L, Moreira D, Talla E, Brochier-Armanet C: A complex cell division 
machinery was present in the last common ancestor of eukaryotes. PLoS 
One 2009, 4(4):e5021. 

38. Cobbe N, Heck MM: The evolution of SMC proteins: phylogenetic analysis 
and structural implications. Mol Biol Evol 2004, 21(2)332-347. 

39. Soppa J: Prokaryotic structural maintenance of chromosomes (SMC) 
proteins: distribution, phylogeny, and comparison with MukBs and 
additional prokaryotic and eukaryotic coiled-coil proteins. Gene 2001, 
278(1-2)253-264 

40. Iyer LM, Burroughs AM, Aravind L: The prokaryotic antecedents of the 
ubiquitin-signaling system and the early evolution of ubiquitin-like beta- 
grasp domains. Genome Biol 2006, 7(7):R60. 

41. De Mot R, Nagy I, Walz J, Baumeister W: Proteasomes and other self- 
compartmentalizing proteases in prokaryotes. Trends Microbiol 1999, 
7(2):88-92. 

42. Darwin KH: Prokaryotic ubiquitin-like protein (Pup), proteasomes and 
pathogenesis. Nat Rev Microbiol 2009, 7(7):485-491. 



Erne ef al. BMC Evolutionary Biology 201 1 , 1 1 :265 Page 1 4 of 1 4 

http://www.biomedcentral.eom/1 471 -2 1 48/1 1 /265 



43. Nunoura T, Takaki Y, Kakuta J, Nishi 5, Sugahara J, Kazama H, Chee GJ, 
Hattori M, Kanai A, Atomi H, ef al: Insights into the evolution of Archaea 
and eukaryotic protein modifier systems revealed by the genome of a 
novel archaeal group. Nucleic Acids Res 2010. 

44. Brochier-Armanet C, Boussau B, Gribaldo S, Forterre P: Mesophilic 
Crenarchaeota: proposal for a third archaeal phylum, the 
Thaumarchaeota. Nat Rev Microbiol 2008, 6(3):245-252. 

45. Bapteste E, Charlebois RL, MacLeod D, Brochier C: The two tempos of 
nuclear pore complex evolution: highly adapting proteins in an ancient 
frozen structure. Genome Biol 2005, 6(1 0):R8S. 

46. Desmond E, Brochier-Armanet C, Forterre P, Gribaldo S: On the last 
common ancestor and early evolution of eukaryotes: reconstructing the 
history of mitochondrial ribosomes. Res Microbiol 2011, 1 62(1 ):53-70. 

47. Desmond E, Gribaldo S: Phylogenomics of sterol synthesis: insights into 
the origin, evolution, and diversity of a key eukaryotic feature. Genome 
Biol Evol 2009,2009:364-381. 

48. Koonin EV: Preview. The incredible expanding ancestor of eukaryotes. 
Ce//2010, 140(5):606-608. 

49. Fritz-Laylin LK, Prochnik SE, Ginger ML, Dacks JB, Carpenter ML, Field MC, 
Kuo A, Paredez A, Chapman J, Pham J, ef al: The genome of Naegleria 
gruberi illuminates early eukaryotic versatility. Cell 2010, 140(5):63 1-642. 

50. Brown JR, Douady CJ, Italia MJ, Marshall WE, Stanhope MJ: Universal trees 
based on large combined protein sequence data sets. Nat Genet 2001, 
28(3):281-285. 

51. Bapteste E, Brinkmann H, Lee JA, Moore DV, Sensen CW, Gordon P, 
Durufle L, Gaasterland T, Lopez P, Muller M, ef al: The analysis of 100 
genes supports the grouping of three highly divergent amoebae: 
Dictyostelium, Entamoeba, and Mastigamoeba. Proc Natl Acad Sci USA 
2002, 99(3):1414-1419. 

52. Brochier C, Forterre P, Gribaldo S: An emerging phylogenetic core of 
Archaea: phytogenies of transcription and translation machineries 
converge following addition of new genome sequences. BMC Evol Biol 

2005, 5(1):36. 

53. Brochier C, Bapteste E, Moreira D, Philippe H: Eubacterial phylogeny based 
on translational apparatus proteins. Trends Genet 2002, 1 8(1 ):1 5. 

54. Rodriguez-Ezpeleta N, Brinkmann H, Burey SC, Roure B, Burger G, 
LofFelhardt W, Bohnert HJ, Philippe H, Lang BF: Monophyly of primary 
photosynthetic eukaryotes: green plants, red algae, and glaucophytes. 
Curr Biol 2005, 1 5(1 4):1 325-1 330. 

55. Burki F, Inagaki Y, Brate J, Archibald JM, Keeling PJ, Cavalier-Smith T, 
Sakaguchi M, Hashimoto T, Horak A, Kumar S, ef al: Large-scale 
phylogenomic analyses reveal that two enigmatic protist lineages, 
telonemia and centroheliozoa, are related to photosynthetic 
chromalveolates. Genome Biol Evol 2009, 1:231-238. 

56. Burki F, Shalchian-Tabrizi K, Pawlowski J: Phylogenomics reveals a new 
'megagroup' including most photosynthetic eukaryotes. Biol Lett 2008, 
4(4):366-369. 

57. Hampl V, Hug L, Leigh JW, Dacks JB, Lang BF, Simpson AG, Roger AJ: 
Phylogenomic analyses support the monophyly of Excavata and resolve 
relationships among eukaryotic "supergroups". Proc Natl Acad Sci USA 
2009, 106(10)3859-3864. 

58. Delsuc F, Brinkmann H, Philippe H: Phylogenomics and the reconstruction 
of the tree of life. Nat Rev Genet 2005, 6(5):361-375. 

59. Rodriguez-Ezpeleta N, Brinkmann H, Roure B, Lartil lot N, Lang BF, 
Philippe H: Detecting and overcoming systematic errors in genome-scale 
phytogenies. Sysf Biol 2007, 56(3):389-399. 

60. Felsenstein J: Cases in which parsimony or compatibility methods will be 
positively misleading. Syst Zool 1978, 27:401-410. 

61. Thomarat F, Vivares CP, Gouy M: Phylogenetic analysis of the complete 
genome sequence of Encephalitozoon cuniculi supports the fungal 
origin of microsporidia and reveals a high frequency of fast-evolving 
genes. J Mol Evol 2004, 59(6)780-791. 

62. Ruiz-Trillo I, Lane CE, Archibald JM, Roger AJ: Insights into the evolutionary 
origin and genome architecture of the unicellular opisthokonts 
Capsaspora owczarzaki and Sphaeroforma arctica. J Eukaryot Microbiol 

2006, 53(5):379-384. 

63. Minge MA, Silberman JD, Orr RJ, Cavalier-Smith T, Shalchian-Tabrizi K, 
Burki F, Skjaeveland A, Jakobsen KS: Evolutionary position of breviate 
amoebae and the primary eukaryote divergence. Proc Biol Sci 2009, 
276(1 657)597-604. 



64. Shalchian-Tabrizi K, Minge MA, Espelund M, Orr R, Ruden T, Jakobsen KS, 
Cavalier-Smith T: Multigene phytogeny of choanozoa and the origin of 
animals. PLoS One 2008, 3(5):e2098. 

65. Ruiz-Trillo I, Roger AJ, Burger G, Gray MW, Lang BF: A phylogenomic 
investigation into the origin of metazoa. Mol Biol Evol 2008, 25(4):664-672. 

66. Baurain D, Brinkmann H, Petersen J, Rodriguez-Ezpeieta N, Stechmann A, 
Demoulin V, Roger AJ, Burger G, Lang BF, Philippe H: Phylogenomic 
evidence for separate acquisition of plastids in cryptophytes, 
haptophytes, and stramenopiles. Mol Biol Evol 2010, 27(7):1 698-1 709. 

67. Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, 
Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein 
database search programs. Nucleic Acids Res 1997, 25(17)3389-3402. 

68. O'Brien EA, Koski LB, Zhang Y, Yang L, Wang E, Gray MW, Burger G, 
Lang BF: TBestDB: a taxonomically broad database of expressed 
sequence tags (ESTs). Nucleic Acids Res 2007, , 35 Database: D445-451. 

69. Katoh K, Misawa K, Kuma K, Miyata T: MAFFT: a novel method for rapid 
multiple sequence alignment based on fast Fourier transform. Nucleic 
Acids Res 2002, 30(14):3059-3066. 

70. Philippe H: MUST, a computer package of Management Utilities for 
Sequences and Trees. Nucleic Acids Res 1993, 21(22)5264-5272. 

71. Durbin R, Eddy SR, Krogh A, Mitchison G: Biological Sequence Analysis. 
Cambridge: University Press Cambridge: 1999. 

72. Finn RD, Mistry J, Schuster-Bockler B, Griffiths-Jones S, Hollich V, 
Lassmann T, Moxon S, Marshall M, Khanna A, Durbin R, ef al: Pfam: clans, 
web tools and services. Nucleic Acids Res 2006, , 34 Database: D247-251. 

73. Jobb G, von Haeseler A, Strimmer K: TREEFINDER: a powerful graphical 
analysis environment for molecular phylogenetics. BMC Evol Biol 2004, 
4:18. 

74. Le SQ, Gascuel O: An improved general amino acid replacement matrix. 
Mol Biol Evol 2008, 25(7):1 307-1 320. 

75. Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference 
under mixed models. Bioinformatics 2003, 1 9(1 2):1 572-1 574. 

76. Lartillot N, Lepage T, Blanquart S: PhyloBayes 3: a Bayesian software 
package for phylogenetic reconstruction and molecular dating. 
Bioinformatics 2009, 25(17)2286-2288. 

77. Lartillot N, Philippe H: A Bayesian mixture model for across-site 
heterogeneities in the amino-acid replacement process. Mol Biol Evol 
2004, 21 (6):1 095-1 109. 

78. Roger AJ, Simpson AG: Evolution: revisiting the root of the eukaryote 
tree. Curr Biol 2009, 1 9(4):R1 65-1 67. 



doi:1 0.1 1 86/1 471 -21 48-1 1 -265 

Cite this article as: Erne ef al: The phylogenomic analysis of the 
anaphase promoting complex and its targets points to complex and 
modern-like control of the cell cycle in the last common ancestor of 
eukaryotes. BMC Evolutionary Biology 201 1 11 :265. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at 
www.biomedcentral.com/submit 



BioMed Central 



